function q = ElnqX_Z(vars)

% E[ln p(X|Z)]

N = size(vars.Z, 1);
D = size(vars.R, 2);
Dlog2pi = D*log(2*pi);

Z = vars.Z;
logdetE = vars.logdetE;

q = -N*D - N*Dlog2pi - sum(sum(Z.*logdetE));
q = q*0.5;
